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Abstract. We develop an approximation theory in Hilbert spaces that gener- 
alizes the classical theory of approximation by entire functions of exponential 
type. The results advance harmonic analysis on manifolds and graphs, thus fa- 
cilitating data representation, compression, denoising and visualization. These 
tasks are of great importance to machine learning, complex data analysis and 
computer vision. 

1. Introduction 

One of the main themes in Analysis is correlation between frequency content 
of a function and its smoothness. In the classical approach the frequency is un- 
derstood in terms of the Fourier transform (or Fourier series) and smoothness is 
described in terms of the Sobolev or Lipshitz and Besov norms. For these notions 
it is well understood [1] , [10] that there exists a perfect balance between the rate of 
approximation by bandlimited functions ( trigonometric polynomials) and smooth- 
ness described by Besov norms. For more recent results of approximations by entire 
functions of exponential type we refer to [5]-[7]- 

The classical concepts and result were generalized to Riemannian manifolds, 
graphs, unitary representations of Lie groups and integral transforms in our work 

ELM, 0- 

The goal of the present article is to develop a form of a Harmonic Analysis 
which holds true in general Hilbert spaces. In Introduction we formulate main 
results obtained in the paper. The exact definitions of all notions are given in the 
text. 

We start with a self-adjoint positive definite operator L in a Hilbert space % and 
consider its positive root D = L 1 / 2 . For the operator D one can introduce notion 
of the Spectral Transform Tu which is an isomorphism between H. and a direct 
integral of Hilbert spaces over R. 

A Paley- Wiener space PW u1 (D),lu > 0, is introduced as the set of all f £ H 
whose image Tof has support in [0,lu]. In the case when H = L2(M.) d and D is 
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a positive square root from the Laplace operator our definition produces regular 
Paley- Wiener spaces of spherical exponential type. 

The domain T> s ,s £ R, of the operator D s ,s £ R, plays the role of the Sobolev 
space and we introduce Besov spaces q = Bj q (D),a > 0, 1 < q < oo, by using 
Peetre's interpolation if -functor [2]. [3], [5], [5], [25]. 

(!•!) B 2, 9 P) = (^^/ 2 )f /r , g) 

where r can be any natural such that < a < r, 1 < g < oo, or < a < r, q = oo. 
It is crucial for us that Besov norms can be described in terms of a modulus of 
continuity constructed in terms of the Schrodinger group e ltD , wave semigroup 
e ltD , or the heat semigroup e~ tD . In what follows the notation || ■ || bellow means 
|| • || fi. We introduce a notion of best approximation 

(1-2) £(./»= M \\f-g\\,fen. 

gePW^iD) 

We also consider the following family of functionals which describe a rate of decay 
of the Spectral transform Tr> 

/ roc \ 1/2 

(1.3) ft(./»=(J \\T D (f)(X)\\ 2 x(x) dm(X)j ,w>0. 

The Plancherel Theorem for To implies that every such functional is exactly the 
best approximation of / by Paley- Wiener functions from PW L0 {D): 

(1-4) K(f,w)=£(f,w)= inf J/-s||. 

Our main results are the following. 

Theorem 1.1. The norm of the Besov space Bj q (D), a>0, l<q<oo is equiv- 
alent to the following norms 

f ds^ 1/q 
(1-5) 11/11+ (/ (s a S(f,s)y- 



(1.6) H/II+ lj2(a ka £(f^ k )) 9 ) ,a>l. 

\k=0 ) 



Jo s 



(1-7) ll/ll 
and 

/ oo \ Vs 

(1.8) II/H+ [Y2(a ka n(f,a k )) q j ,a>l. 

Theorem 1.2. A vector f £ H belongs to Bf q (D),a > 0, 1 < g < oo, i/ and onfo/ 
z/ i/iere exists a sequence of vectors fk = fk{f) € PW a k {D), a > 1, fc G N smc/i £/ia£ 
series fk converges to f in H and the following inequalities hold for some 
c\ > 0,C2 > which are independent on f £ BJ^-D) 

/ oo \ 1 /l 

(1-9) ciII/IIbj,,^) < £ ( afe "IIMI) 9 < C8ll/llB ?i ,(u),a > 1. 

\fe=0 / 
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In the case when a > 0, q = oo one has to make appropriate modifications in the 
above formulas. 

According to (11.41) the functional £(f,ui) is a measure of decay of the Spectral 
Transform To and the Theorems 11.11 and 11.21 show that Besov spaces on a manifold 
M describe decay of the Spectral transform J-jj associated with any appropriate 
operator D. 

In the case Ti = L2(M d ) the Theorems 11.11 and 1 1 . 21 are classical and can be found 
in PQ, [TU] and [32]. In the case when T-L is L2-space on a Riemannian manifold 
or a graph and D is the square root from the corresponding Laplace operator the 
Theorems 1 and 2 were proved in our papers [Tl]-[20]. 

2. Paley- Wiener subspaces generated by a self-adjoint operator in a 

HlLBERT SPACE 

Now we describe Paley- Wiener functions for a self-adjoint positive definite oper- 
ator D in Ti. According to the spectral theory [3] for any self-adjoint operator D in 
a Hilbert space T-L there exist a direct integral of Hilbert spaces X = J X(X)dm(X) 
and a unitary operator Td from T-L onto X, which transforms domain of D k , k 6 N, 
onto Xk = {x 6 X\\ k x 6 X} with norm 

/ poo \ 1/2 

(2.1) \\x(X)\\ Xk = U \ 2k \\x(X)\\l w dm(X)j 

besides J r n(D k /) = \ k (J~Df), if / belongs to the domain of D k . As it is known, X 
is the set of all m-measurable functions A — > x(X) € X(\), for which the norm 

/ i-oo \ 1/2 

Mix = (j o \\x(X)f xW dm(X)) 

is finite. 

Definition 1. We will say that a vector / from T-L belongs to the Paley- Wiener 
space PW U1 (D) if the support of the Spectral transform Fuf belong to [0,uj]. For 
a vector / g PW U (D) the notation uif will be used for a positive number such that 
[0,u;/] is the smallest interval which contains the support of the Spectral transform 

Using the spectral resolution of identity P\ we define the unitary group of oper- 
ators by the formula 

/■oo 

e ltD f= / e u ^dP T f,f €H,t€R. 
Jo 

Let us introduce the operator 

(2-2) RS/ = W^W^ (k ' im)DfJ € n ^ >0 - 

Since ||e lt£ /|| = ||/|| and 



(2.3) 
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the series in (|2.2[) is convergent and it shows that R£, is a bounded operator in % 
with the norm uj: 

(2.4) ||B£/||<«ll/ll,/€tt. 

The next theorem contains generalizations of several results from the classical 
harmonic analysis (in particular the Paley- Wiener theorem) and it follows essen- 
tially from our results in [13], [14], [19] . 

Theorem 2.1. The following statements hold: 

(1) the set LL>o PWui(D) is dense in %; 

(2) the space PW L0 (D) is a linear closed subspace in W; 

(3) a function f £ % belongs to PW bJ {D) if and only if it belongs to the set 

oo 

2?oo = f) T> k (D), 

k=l 

and for all s £ M + the following Bernstein inequality takes place 
(2-5) \\D s f\\<uj s \\f\\; 

(4) a vector f £ H belongs to the space PW Uf {D), < ujf < oo, if and only if 
f belongs to the set T>oc, the limit 



lim \\D k f\\ 1/k 

exists and 



(2.6) lim \\D k f\\ 1 / k =u f . 

(5) a vector f £ H belongs to PW UJ (D) if and only if f £ T)^ and the upper 
bound 

(2.7) sup {uj- k \\D k f\\) < oo 

keN 

is finite, 

(6) a vector f € % belongs to PW U {D) if and only if f € and 

(2-8) lim^ 00 ||^/H 1/fc =^<^. 

In this case uj = ujj. 

(7) a vector f £ H belongs to PWuj(D) if and only if it belongs to the to the 
set 2?oo and the following Riesz interpolation formula holds 

(2.9) (iD) n f = (B%) n f,neN; 

(8) / € PW U] {D) if and only if for every g 6 % the scalar-valued function of 
the real variable (e ltD f,g} ,t £ R 1 , is bounded on the real line and has an 
extension to the complex plane as an entire function of the exponential type 
w; 

(9) / G PWu(D) if and only if the abstract-valued function e ltD f is bounded 
on the real line and has an extension to the complex plane as an entire 
function of the exponential type u); 

(10) / £ PWcj(D) if and only if the solution u(t), t £ M 1 of the Cauchy problem 
for the corresponding abstract Schrodinger equation 

.du{t) 



t)j =Du(t),u(0) = f,i = V = l ) 
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has analytic extension u(z) to the complex plane C as an entire function 
and satisfies the estimate 

\Hz)\\ H <e^\\f\\u. 

3. Direct and Inverse Approximation Theorems 

Now we are going to use the notion of the best approximation (|1.2p to introduce 
Approximation spaces E^ g (Z?),0 < a < r,r £ N.l < q < oo, as spaces for which 
the following norm is finite 

(3-1) = 11/11 + QT (s a £(f,s)) q ^ 

where < a < r,l < q < oo, or < a < r, q — oo. It is easy to verify that this 
norm is equivalent to the following "discrete" norm 

(3-2) II/H+ (]T( a ^(/V)n ,a>l, 

The Plancherel Theorem for J~d also gives the following inequality 

/ /-oo \ 1/2 

(3.3) £(f,u)<u- k i^J \\T D (D k f)(\)\\ 2 x{x) dm(\)) < u;- k \\D k f\\. 

In the classical Approximation theory the Direct and Inverse Theorems give 
equivalence of the Approximation and Besov spaces. Our goal is to extend these 
results to a more general setting. 

For any / £ T-L we introduce a difference operator of order m £ N as 

m 

(3.4) A™/ = (-l) m+1 ^2(-iy- 1 Cl l e^ iD ^f,T£ R. 

3=0 

and the modulus of continuity is defined as 

(3-5) n m (f,s)= sup ||A-/|| 

|t|<s 

The following Theorem is a generalization of the classical Direct Approximation 
Theorem by entire functions of exponential type flOl . 

Theorem 3.1. There exists a constant C > such that for all co > and all f 

(3.6) £(f, u) < ^rQm-k (D k f, 1/w) , < k < m. 
In particular the following embeddings hold true 

(3.7) Bl q (D)cE«(D),l<q<oo. 

Proof. If h £ Li(M) is an entire function of exponential type uj then for any f £% 
the vector 

/oo 
h(t)e ltD fdt 
-oo 

belongs to PW U (D). Indeed, for every real r we have 



e 



irD g= / h(t)e t(t+T)D fdt= / h(t-T)e ltD fdt 
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Using this formula we can extend the abstract function e lTD g to the complex plane 
as 

/oo 
h(t - z)e itD fdt. 
-OO 

Since by assumption h 6 Li(R) is an entire function of exponential type oj we have 

/>oo pOO 
\h(t-z)\dt<\\f\\e?\'\ / \h{t)\dt. 
-oo J —oo 

It shows that for every functional g* G H the function {e LzD g,g*) is an entire 
function and 



/oo 
\h{t)\dt. 
-oo 



In other words (e lzD g,g*^ is an entire function of the exponential type uj which is 
bounded on the real line and application of the classical Bernstein theorem gives 
the following inequality 



(I)' <«"«•»•> 



<^sup\(e" D g,g-)\. 



Since 



we obtain for t = 



/OO 
\h(r)\dr. 
-oo 

Choosing g* such that ||^*|| = 1 and (D k g,g*) = \\D k g\\ we obtain the following 
inequality 

/oo 
\h(r)\dr 
-oo 

which implies that g belongs to PW U1 {D). 
Let 



(3.8) h(t) = a 

where n > m + 3 is an even integer and 



sin(i/n)\™ 



a = 



sin(t/n) 



dt 



With such choice of a and n the function ft will have the following properties: 

(1) h is an even nonnegative entire function of exponential type one; 

(2) h belongs to and its Li(R)-norm is 1; 

(3) the integral 



(3.9) 



/OO 
h(t)\ 
-oo 



t\ m dt 



is finite. 
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Consider the following vector 

(3.10) Qr n (/) = r ^{(-lr^A^z+z}^ 



where 

m m 

(3.11) (-l) m+1 A™/ = (-l) m+1 ^(-l) J '- 1 C4e is ( iD )/ ^/y ' 8 ( iD / - /, 

3=0 j=i 

and 

(3.12) & x + 6 2 + ... + 6 m = l. 
The formulas (|3.10[) and p. Ill) imply the following formula 

/oo m />oo 

h(t)J2b j e ji(W) fdt= / $We t(lD) M. 
- oo j — 1 J —oo 

where 

.<.)-£>(?)*(.? 

Since the function /i(t) is of the exponential type one every function h(tu>/j) is of 
the type w/j. It also shows that the function $(£) is of the exponential type ui as 
well. 

Now we estimate the error of approximation of / by Q h ' (/). If the modulus of 
continuity is defined as 

(3.13) n ro (/,*)= sup |K7|| 

\T\<8 

then since by (|3 . 1 0|) 

/oo 
-OO 

we obtain 

/OO pOO 
h(t) AT/ J dt< / h(t)n m (f,t/u)dt. 
-oo J —oo 

Now we are going to use the following inequalities 

(3.14) n m (f, s ) < s k n rn _ k (D k f,s) 



(3.15) a m (/, as) < (1 + a) m fi m (/, s), a G M+. 
The first one follows from the identity 

(3.16) A k f = (e UD -lff= [ ... / e^ + -^ D D k fdr 1 ...dr k , 

Jo Jo 

where / is the identity operator and k 6 N. The second one follows from the 
property 

fii(/,«i + *2)<ni(/,si) + fii(/,« 2 ) 

which is easy to verify. We can continue our estimation of E(f,uj). 

E{f,uj)< h(t)n m (f,t/oj)dt< m ~ kK k 1 > h(t)\t\ k (i + \t\y n - k dt< 



BESOV VECTORS AND PALEY- WIENER VECTORS IN HILBERT SPACES 



UJ 



-n m - k (D k f,i/u), 



where the integral 

cLk= I Ht)\t\ k {i + \t\r~ k dt 



is finite by the choice of h. The inequality (|3 . 6[) is proved and it implies the second 
part of the Theorem. □ 

In fact we proved a little bit more. Namely for the same choice of the function 
h the following holds. 

Corollary 3.1. For any < k < m, k, m G N, here exists a constant CiL k such 
that for all < U) < oo and all f G % the following inequality holds 

(3.17) £{f, w) < || Qt m {f) - /|| < %^O ro _ fc (D h f, l/ w ) , 

/oo 
h(*)|t|*(l + |t|rdt,0<fe<m ) 
-oo 

and i/ie operator 

Q% m :H->PW U (D) 

is defined in i3.10\) . 

Next, we are going to obtain the Inverse Approximation Theorem in the case 
q = oo. 

Lemma 1. If there exist r > a — n > 0, a > 0, r, n G N, such that the quantity 

(3.18) b« (/) = sup ( S "~ q O. (£>"/, «)) 

is finite, then there exists a constant A = A(n, r) for which 

(3.19) S ups a 8(f,s)<A(n,r)b^ ntr (f)- 

s>0 

Proof. Assume that (|3.18[) holds, then 

n r (D n f, s )<b^ r (f)s a - n 

and (|3.17l) implies 

(3.20) A(n,r)b^„, r (./>- Q . 

Lemma is proved. □ 
Lemma 2. If for an f G % and /or an a > the following upper bound is finite 

(3.21) sup S a £(/, s ) =T(/,a) < oo, 

s>0 

i/ien /or every r > a — n > 0, a > 0, r, n G N, i/iere exists a constant C{a, n, r) suc/i 
i/iai i/ie next inequality holds 

(3-22) b» ,„-(/) < C(a,n,r) (||/|| +T(/,a)) . 
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Proof. The assumption implies that for a given / £ H and a sequence of numbers 
a J , a > 1, j = 0, 1, 2, ... one can find a sequence c/j £ PW a j (D) such that 

||/-<fe||<T(/ J a)a-'' Q ,a>l. 

Then for 

(3-23) f = go, U = 9j - 9j-i G PW 8 i (£>), 

the series 

(3-24) / = /o + /i+/ 2 + .... 

converges in T~L. Moreover, we have the following estimates 

ll/oll =W< 11.90- /II + ll/ll < ll/ll +T(/,a), 

11/,'H < \\f-9i\\ + \\f-9j-i\\ <T(f,a)a-i a +T(f,a)a-U-V a =T(f,a)(l+a a )a-i a , 
which imply the following inequality 

(3-25) \\f j \\<C(a,a)a-^(\\f\\+T(f,a)),j£N. 

Since fj £ PW a j (D) we have for any n £ N 

(3-26) \\D n fj\\<a^\\fj\\,a>l, 

we obtain 

\\D n fi\\ < C(a,a)a-* a - n ) (||/|| +T(f,a)) 
which shows that the series 

je® 

converges in % and because the operator D n is closed the sum / of this series 
belongs to the domain of D n and 

D n f = Y / D n f 3 . 

jew 

Next, let Fj = D n fj then we have that D n f = £\ Fj, where Fj € PW a ,(D) and 
according to (f3~25]l and (pT26|) 

(3.27) ll^ll = H^/,-11 < a^\\fj\\ < C(a,a)a-^- n H\\f\\+T(f,a)). 
Pick a positive t and a natural N such that 

(3.28) a~ N <t <a- N+1 ,a> 1, 

then we obviously have the following formula for any natural r 

N-l oo 

(3.29) A r t D n f = & r t Fj + £ & r tFj> 

j=0 j=N 

where A[ is defined in p.4[) . Note, that the Bernstein inequality and the formula 
(l3~16l) imply that if / £ PW U {D), then 

(3-30) UAf/H < 

Since (|3.27p and p.28j) hold we obtain for j < N — 1 the following inequalities 

IIA^-H < (a?ty\\F s \\ < C(a,a)(\\f\\+T(f,a))a^ +r - a ^ N -^\a> 1. 
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These inequalities imply 
w-i 



J2 A t F i 



3=0 



<C(a,a)(\\f\\+T(f,a))a 



N-l 

-r(N-l) a (n+r-a)j _ 

3=0 



1 _ a {n+r- a )N 



C(a,a)(||/||+T(/,a))a-^-^ T - < 

(3-31) C(a,a,n,r)(\\f\\+T(f,a))t a - n . 
By applying the following inequality 

||A^-||< 21^-11 
to terms with j > N we can continue our estimation as 



J2 A t F 3 



j=N 



<2 r C(a,a)(\\f\\+T(f,a))J2a- (a - nh 

j=N 



C(a,a)2 r (\\f\\+T(f,a))a- N ^- n \l-a (n -^)- 1 < 

(3-32) C(a,a,n,r)(\\f\\+T(f,a))t a - n . 
It gives the following inequality 

\\A r t D n f\\ < C(a,a,n,r)t a - n (\\f\\+T(f,a)), 
from which one has 

O- (D n f, s) < C(a, a, n, r)(||/|| + T(/, a))s a ~ n , s > 0, 

and 

b^ n!r (f)<C(a,a,n,r)(\\f\\+T(f,a)). 
The Lemma is proved. 

Our main result concerning spaces B2 ^(-D^a > 0, is the following. 



□ 



Theorem 3.2. The norm of the space B^D), a > 0, is equivalent to the following 
norms 



(3.33) 
(3.34) 
(3.35) 
(3.36) 



\f\\+sup(s a £(f,s)), 
f\\+ sup (s a 1l(f,s))), 

s>0 

I- sup (a ka £(f,a k )) ,a>l, 

k&i 

sup (a ko "Jl(f, a k ))) , a > 1. 



Moreover, a vector f <E Ti belongs to BjocX-D), a > 0, if and only if there exists 
a sequence of vectors = fk(f) € /W a jc(D),a > 1, such that the series fk 
converges to f in % and 

(3.37) 



ci||/||b ?oo( z>) <sup(a to ||M|) <c 2 ||/||B ?0O (2 J ),a>l, 

fceN 



for certain c\ = ci(D,a),c 2 = c 2 (D 7 a) which are independent of f € B^^-D). 
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Proof. That the norm of oo(D), a > 0, is equivalent to any of the norms f|3 . 33(> - 
(|3.36p follows from the last two Lemmas and (|1.4p . 

Next, if the norm (13.33)) is finite then it was shown in the proof of the last Lemma 
that there exists a sequence of vectors fk — fk(f) G PW a k(D), a > 1, such that the 
series fk converges to / in H. Moreover, the inequality (|3.25p shows existence of 
constant c which is independent of / € Ti?! oo(D) for which the following inequality 
holds 

sup(a fcQ ||M|) <c||/|| B - (D),a> 1, 
fceN 

Conversely, let us assume that there exists a sequence of vectors fk = fk(f) € 
PW a k(D), a > 1, such that the series Yl fk converges to / in TL and 

sup(a ta ||/ fe ||) <oc. 



We have 



£{f,a 



N \ 



< 



N-l 

£. 

k=0 



£iia 



k=N 



C sup {a ka \\f k \ 



< sup (a 

keN 
\ „-Na 



IMI) X>-°°*< 



k=N 



Since we also have 



supa JVQ £(/,a JV ) < Csup (a fcQ ||/ fe | 

N keN 



/II <EHMI <sup( a fcQ ||/ fc ||)^a- Qfc , a >l, 



k-e 



the Theorem is proved. 



□ 



The Theorems 11.11 and 11.21 from Introduction are extensions of the Theorem 13.2 
to all indices 1 < q < oo. Their proofs go essentially along the same lines as the 
proof of the last Theorem and are omitted. 
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